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PREDICTING  WHEELED  VEHICLE  MOTION 
RESISTANCE  IN  SHALLOW  SNOW 

George  Blaisdell 


INTRODUCTION 

A  vehicle  traveling  through  snow  is  required  to  expend  a  greater 
amount  of  energy  than  it  would  on  an  undeformable  surface*  Since  this 
energy  difference  is  visually  accounted  for  by  the  formation  of  a  rut, 
many  researchers  have  postulated  that  the  energy  required  to  compact  a  snow 
layer  should  be  equal  to  the  increase  in  motion  resistance  of  a  vehicle 
moving  in  that  snow  layer.  Provided  the  energy  density  of  a  snow  layer  can 
be  determined,  the  compaction  energy,  and  hence  the  motion  resistance,  can 
be  calculated  for  a  wheel  or  track  moving  through  snow. 

PAST  WORK 

The  first  attempt  to  use  an  energetics  approach  for  the  prediction 
of  vehicle  mobility  over  snow  was  made  by  Landauer  and  Royse  (1956),  who 
generated  force  vs  sinkage  curves  for  a  natural  snow  with  a  hydraulically 
driven  plate.  By  integrating  the  curves  using  planimetry,  they  established 
an  energy  per  unit  area  term,  w.  When  compared  with  motion  resistance 
values  measured  in  the  same  area  with  a  towed  tire,  the  compaction  energy 
determined  by  their  method  yielded  lower  values.  The  authors  thus 
concluded  that  the  cause  of  increased  motion  resistance  had  little  to  do 
with  the  work  of  compaction. 

Bekker  (1956)  and  Nuttall  (1957)  continued  the  idea  of  equating 
motion  resistance  to  the  work  of  compaction.  Each  established  a  mathe¬ 
matical  relationship  among  pressure,  density  and  sinkage.  Upon  evaluation 
of  the  Bernstein  (1913)  equation 

z 

R  =  b  /  °p  dz  (1) 

o 


where  R  =  motion  resistance 
b  =  contact  width 
z0  =  sinkage 
p  ®  pressure, 

Bekker  (1956)  and  Nuttall  (1957),  respectively,  express  motion  resistance  as 


R  = 


n  +  1 


(2) 


V 


Figure  1.  Hyperbolic  pressure 
distribution  along  the  wheel/snow 
interface  (Liston  1974) . 


where  k  is  a  coefficient  of  proportionality  and  n  a  function  of  snow 
properties,  and 

hi  P  Y. 

R  =  F  t  r  (1  ‘  t*  (3) 

*  o  Yf 

where  F  =  applied  vertical  force 

hj  =  initial  snow  depth 
£  =  contact  length 
Pz  =  pressure  at  sinkage  z 
PD  *  average  ground  pressure 
■Yi  =  initial  snow  density 
Yf  =  final  snow  density. 

In  a  slightly  different  approach,  Liston  (1974)  defined  the  work  E 
expended  in  compacting  snow  from  its  initial  volume  Vj  to  a  final  volume 
V2  as 

v2 

E  =  /  p  dV,  (4) 

▼l 

and  assumed  that  the  pressure  vs  volume  relationship  varies  hyperboli- 
cally  (Fig.  1).  This  relationship  is  based  on  an  analysis  of  the  wheel/ 
snow  interface  and  assumes  that  1)  the  only  dimensional  change  that  occurs 
during  compaction  from  Vi  to  V2  is  in  the  snow  depth  and  2)  the  pressure 
along  the  wheel/snow  interface  equals  the  atmospheric  pressure  at  the  point 
of  contact  between  the  snow  surface  and  the  wheel  and  equals  the  inflation 
pressure  at  the  lowest  point  of  contact.  The  pressure-volume  relationship 
then  becomes 


P  " 


(5) 


where  p^  is  tire  inflation  pressure.  Introducing  eq  5  into  eq  4  and 
completing  the  integration  yields 


2 


(6) 


E  -  Pl»2  1“  ^  ■ 

By  defining  the  force  resisting  motion  as  R,  its  relation  to  compaction  is 
then 


E  -  RA  • 


(7) 


Noting  that  the  volume  V  is  the  product  of  contact  length  A,  contact  width 
b  and  the  snow  depth  h,  and  combining  eq  6  and  7,  the  resistance  equation 
becomes 


R  =  Pibh2  In 


(8) 


where  hi  is  initial  snow  depth  and  h£  is  final  snow  depth  after  vehicle 
passage. 

The  Bekker  (1956) ,  Nuttall  (1957)  and  Liston  (1974)  approaches  are  all 
similar  in  that  they  equate  the  work  of  compaction  of  snow  to  the  product 
of  distance  traveled  and  force  resisting  motion.  Likewise,  each  attempts 
to  develop  one  mathematical  expression  that  describes  the  force  vs  sinkage 
relationship  for  all  snows.  Realizing  that  the  effects  of  temperature, 
density,  grain  size,  stratification,  etc.  combine  to  produce  different 
force  vs  sinkage  curves  in  different  snows,  Harrison  (1975)  suggests 
obtaining  the  energy  of  compaction  from  an  actual  force  vs  sinkage  record. 
All  previous  approaches  require  information  obtained  in  the  field  anyway, 
so  this  does  not  pose  any  new  inconvenience.  He  then  suggests  that  an 
energy  per  unit  volume  term  be  used  to  characterize  the  snow  to  calculate 
motion  resistance. 

Harrison  notes  that  the  deformation  z  nears  a  maximum  value  ze  as 
maximum  compaction  is  approached  during  a  compression  test  with  snow  con¬ 
fined  on  top  and  bottom  only.  This  is  considered  a  shallow  snow  case 
because  the  pressure  bulb  extends  to  the  bearing  surface  under  the  snow. 
Further  loading  of  the  snow  once  critical  density  yCr  is  reached 
(occurring  at  sinkage  ze)  causes  plastic  flow  (mass  movement  of  the  snow 
rather  than  further  compaction)  to  take  place.  It  has  been  found  that  a 
critical  density  does  exist  for  snow  and  that  the  majority  of  vehicles 
traveling  through  shallow  snow  do  compact  the  snow  to  this  critical  density 
(Bennett  1973,  Hinic  1974,  Kinoshita  and  Akitaya  1970). 
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The  area  under  the  force  vs  sinkage  curve  between  the  limits  of  zero 
sinkage  and  ze  is,  then,  equal  to  the  total  work  E  expended  in  compacting 
the  snow  to  its  critical  density.  To  obtain  an  energy  per  unit  volume  term 
to  the  total  work  E  is  divided  by  the  total  volume  prior  to  loading  that 
became  involved  in  compaction: 


This  equation  assumes  that  all  the  snow  involved  in  compaction  has  been 
brought  to  its  critical  density  and  that  this  energy  term  is  independent  of 
the  volume  used  to  obtain  it.  Variations  from  these  assumptions  and  their 
possible  effects  will  be  discussed  later. 

Harrison  (1975,  1981)  has  suggested  the  use  of  a  circular  plate- 
sinkage  device  to  obtain  values  of  w.  Because  of  its  simple  geometry  and 
the  nature  of  the  pressure  bulb  developed  beneath  the  plate,  calculation  of 
the  volume  of  material  undergoing  compaction  is  easily  accomplished. 

Figure  2  shows  both  thin-section  and  snow  pit  views  of  the  pressure  bulb 
beneath  a  plate.  It  is  readily  apparent  that  the  volume  of  snow  disturbed 
is  confined  to  a  cylinder  with  the  diameter  of  the  plate  used.  Since  we 
are  concerned  with  a  shallow  snow,  the  cylinder's  height  is  the  snow 
depth.  Thus,  the  expression  for  energy  of  compaction  for  a  plate-sinkage 
device  is 


a) 


E 


A  hi 
P  1 


(10) 


where  Ap  is  the  area  of  the  plate  and  Vq  is  initial  snow  depth. 

The  E  term  in  eq  9  and  10  is  found  from  a  plot  of  force  vs  sinkage 
obtained  during  the  plate  test  (Fig.  3).  Using  planimetry,  the  area  under 
the  curve  is  found  between  the  limits  of  z  ■  0  and  z  ■  ze.  The  quantity 
ze  may  be  found  from  the  force-sinkage  curve  where  the  curve  becomes 
asymptotic  or  by  calculation  from  (Harrison  1973) 


z  -  hi  (1  -  — )  (ID 

e  1  Y 

1  cr 

where  Yi  is  initial  snow  density  and  Ycr  =  critical  density. 

It  should  be  cautioned  that  the  pressure  bulb  must  intersect  the 
ground  surface  in  order  to  use  h^  in  eq  10  as  the  pre- compaction  height  of 
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Figure  3.  Idealized  force  vs  sinkage 
curve  generated  by  a  plate-sinkage 
device  in  shallow  snow. 


material  involved  in  densif ication.  Harrison  (1981)  also  recommends  a 
minimum  plate  diameter  of  20  cm  be  used  to  minimize  size  effects. 

Calculation  of  u)  using  other  devices  is  not  precluded  but  may 
introduce  complications  in  the  determination  of  both  E  and  V.  Unlike  the 
plate  test,  other  tests  may  not  compact  all  of  the  material  disturbed  by 
the  loading  device  to  its  critical  density.  The  volume  disturbed  by  other 
pressure-sinkage  devices  may  be  of  an  irregular  shape  and  thus  difficult  to 
compute.  In  addition,  one  must  decide  the  percentage  of  the  volume  not  at 
Ycr  to  add  to  that  at  ycr  for  the  volume  term  in  eq  9.  In  determining 
E,  one  must  choose  a  point  along  the  load  vs  deformation  curve  (generated 
by  the  device)  that  corresponds  with  having  reached  critical  density.  This 
point  should  also  represent  the  onset  of  plastic  flow  in  the  snow. 

To  obtain  the  resistance  to  motion  R  from  snow  compaction  for  one  tire 
of  a  vehicle  traveling  through  shallow  snow,  Harrison  (1975)  suggests  that 
the  energy  per  unit  volume  term  be  multiplied  by  the  volume  of  snow  which 
the  vehicle  compacts,  or 

R  =  bhiu  (12) 

where  b  is  the  track  width  or  effective  tire  width  and  to  the  energy 
required  per  unit  volume  to  compact  snow  to  its  critical  density. 

An  approach  similar  to  Harrison's  was  used  by  Selig  (1971)  for 
evaluating  soil  compaction  equipment.  For  a  towed  rigid  wheel  roller, 
the  compactive  effort  e  is  defined  as 
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(13) 


e  -  FtLP 

where  F^-  =  average  towing  force  (equivalent  to  motion  resistance 
in  the  previous  approaches) 

L  =  distance  towed 

P  =  number  of  passes  over  distance  L. 

This  relationship  indicates  that  the  compact ive  effort  is  provided  by 
the  towing  unit  through  the  drawbar.  Selig  suggests,  however,  that  the 
towing  force  can  be  related  to  the  roller  weight  by 

F  =  fW  (14) 


where  f  is  the  coefficient  of  rolling  resistance  and  W  the  total  roller 
weight. 

The  volume  of  soil  Vc  considered  to  be  absorbing  the  compactive 
effort  is  defined  as 

V  =  LBt  (15) 

c 


where  B  is  the  roller  width  and  t  the  layer  thickness  after  compaction. 

Selig  uses  the  compacted  thickness  rather  than  loose  thickness  to  be 
consistent  with  the  convention  used  in  the  standard  laboratory  impact 
compaction  test. 

Defining  a  compactive  effort  per  unit  volume  term  Ec  by  combining 
eq  13-15  yields 


E 

c 


fWP 

Bt 


(16) 


For  pneumatic  tire  rollers  this  relationship  basically  remains  the  same  ex¬ 
cept  that  B  is  taken  as  the  sum  of  all  the  tire  widths  (provided  the  tire 
separation  is  greater  than  the  tire  width) . 

Rearranging  eq  16  to  predict  average  towing  force  Ft  we  see  that 

F  =  fW  =  E  Bt  (17) 

t  c 

for  one  pass  of  the  vehicle  (P  =  1).  This  relationship  closely  resembles 
Harrison's  expression  for  motion  resistance  (eq  12).  The  primary 
^■1^®*-®^®®  between  eq  17  and  12  lies  in  the  height  of  material  used  in  the 
volume  calculation,  one  using  pre-compaction  height  (Harrison)  and  the 
other  using  the  post-compaction  height  (Selig) . 
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Ec  and  u>  are  clearly  similar  energy  terms  and  are  found  in  a  like 
manner.  As  stated  previously,  w  is  found  from  a  plate-sinkage  test  where  a 
flat  circular  plate  is  forced  into  the  snow  along  a  straight  path  roughly 
perpendicular  to  the  snow  and  ground  surfaces.  The  value  of  Ec  is  based 
on  the  relationship  of  the  measured  soil  densities  in  the  field  to 
densities  produced  by  the  standard  and  modified  laboratory  compaction 
tests  (Wu  1976).  Both  Ec  and  m  are,  then,  determined  from  the  application 
of  a  force  perpendicular  to  the  ground  surface.  Both  approaches,  though, 
use  this  energy  term  to  predict  a  horizontal  motion  resisting  force. 

Selig  attempts  through  f  in  eq  14  to  relate  the  force  required  to 
advance  the  wheel  with  the  vertical  load  acting  on  that  wheel.  However,  f 
is  an  experimentally  determined  value  and  thus  requires  running  tests  in 
which  the  motion  resistance  Ft  is  measured.  From  these  tests,  tables  of 
f  values  are  generated  for  various  types  of  soil  and  compaction  equipment. 
Therefore,  f  is  an  empirically  derived  correction  factor  which  relates 
motion  resistance  to  soil  and  vehicle  characteristics.  Attempts  to  predict 
motion  resistance  using  eq  17  require  relying  on  the  completeness  and 
applicability  of  tabulated  values  of  f. 

Using  Harrison’s  approach  for  making  predictions  of  motion  resistance, 
a  plate-sinkage  test  may  be  performed  in  the  snow  of  interest  from  which  an 
0)  value  is  obtained.  The  (0  value  can  then  be  used  with  eq  12  to  transfer 
the  information  gained  from  a  vertical  application  of  force  to  a  horizontal 
motion  resisting  force. 

THE  MECHANICS  OF  A  WHEEL  TRAVELING  THROUGH  SNOW 

Each  of  the  six  previous  approaches  to  predicting  motion  resistance 
involves  calculations  based  on  a  vertically  applied  force.  Yet,  motion 
resistance,  as  a  vector  quantity,  is  directed  horizontally.  Both  Harrison 
(1981)  and  Selig  (1971)  recognize  this  fact  and  attempt  to  account  for 

it  through  an  empirical  approach.  It  should  be  possible,  however,  to 
mathematically  describe  the  magnitude  and  direction  of  the  forces  at  the 
tire/material  interface. 

The  wheel  can  be  assumed  to  be  constructed  of  an  infinite  number  of 
flat  plates  arranged  tangentially  to  the  circumference.  As  the  wheel  moves 
forward  each  plate  is  pushed  into  the  snow  along  a  path  described  by  the 
motion  of  a  point  on  the  tire  surface  (Fig.  4a).  A  pressure  bulb  develops 
beneath  each  flat  section  of  the  tire  which  grows  as  distance  along  the 
sinkage  path  increases  from  the  initial  tire/snow  contact  point  (Fig.  4b). 
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b. 

Figure  4.  Sinkage  path  (a)  and  pressure  bulb 
development  (b)  beneath  a  wheel  made  of  plates. 

The  pressure  bulb  continues  to  grow  until  the  influence  of  the  ground  or 
bearing  surface  is  felt,  at  which  point  the  snow  becomes  compressed  to  its 
critical  density  and  no  further  sinkage  takes  place •  The  volume  of  snow 
compacted  by  a  wheel  in  moving  from  position  A  to  position  B  (Fig.  4b)  can 
be  found  from 

V  =  m  b  [s  +  (h  -  z)]  (18) 

where  m  is  the  arc  length  along  the  tire/snow  contact  surface  and  s  the 
length  of  the  sinkage  path. 

The  energy  required  to  move  the  wheel  from  position  A  to  B  is  the 
product  of  the  force  necessary  to  advance  the  compacting  tire  and  the  dis¬ 
tance  traveled.  This  force  has  two  components.  One  component  is  developed 
parallel  to  the  snow/vehicle  contact  surface  by  the  driven  wheel.  The 
other  component  is  the  portion  of  the  vehicle  weight  acting  on  the  wheel  of 
interest  and  directed  vertically. 

The  forces  resisting  forward  motion  can  be  found  by  going  back  to 
the  wheel  made  up  of  plates.  If  one  assumes  that  the  force  and  deformation 
vectors  coincide  as  the  snow  is  compacted,  each  plate  feels  a  resistance  to 
sinkage  FTi  which  acts  tangentially  to  the  sinkage  path  at  the  plate/ snow 
contact  point  (Fig.  5).  The  magnitude  of  this  resistance  is  a  function  of 
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Critical  Density 


Figure  5.  Force  vectors  resisting  forward  wheel 
motion  in  snow. 


the  plate’s  sinkage  along  the  curved  path,  just  as  in  the  case  of  a  simple 
plate-sinkage  test.  The  vector  sum  of  these  forces  (acting  on  all  the 
plates)  represents  the  total  resistance  to  compaction  FT  offered  by  the 
snow.  The  horizontal  component  F^  of  this  resultant  force  represents 
the  true  resistance  to  forward  motion  that  is  opposing  the  energy  input  by 
the  driving  wheels.  If  one  ignores  traction,  it  is  only  this  horizontal 
component  of  force  that  makes  travel  over  deformable  surfaces  consume  more 
energy  than  movement  over  non-deforming  materials. 


CALCULATION  OF  THE  FORCE  NECESSARY  TO  COMPACT  SNOW  WITH  A  WHEEL 

To  determine  the  magnitude  and  direction  of  the  resultant  force  Fj, 
the  force  vs  sinkage  relationship  along  a  curved  sinkage  path  and  the 
length  of  the  sinkage  path  is  needed.  This  curved  path,  called  a  tro¬ 
choid,  is  generated  by  describing  the  motion  of  a  point  on  the  circum¬ 
ference  of  the  wheel  as  it  rotates.  The  general  form  of  the  equation  for  a 
trochoid  describing  wheel  motion  is 

x  =  a0  -  r  sine 

y  =  r(l  -  cos6)  (19) 

where  r  =  wheel  radius 

a  =  radius  of  a  hypothetical  wheel  which  would  cover  the  same 
distance  as  the  real  wheel  after  moving  through 
angle  0,  but  without  any  slippage  (a  =  r  for  a  wheel  with  no 
slippage) 

0  =  angle  (in  radians)  between  a  vertical  radius  and  the 

radius  which  ends  at  the  point  under  consideration. 
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The  angle  @a  for  the  point  of  initial  tire  contact  with  the  snow  surface 
(Fig.  4a)  can  be  found  from  sinkage  by 

r 1  -z 

0A  =  2tt  -  arccos  — —  (20) 

where  r1  is  the  effective  radius  (rf  «  r  -  6  where  6  is  tire  deflection). 
If  we  begin  by  assuming  a  non-driven  wheel  with  no  slip,  a  =  r  in  eq  19. 
The  length  of  the  sinkage  path  can  then  be  found  from 

2tt 

A 


s  =  / 

niQv  * 


which  yields 

2tt 

s  =  -4r  cos  7T  I  .  (22) 

max  2  I 

0A 


It  is  along  this  distance  s  that  the  force  required  to  compact  snow 
varies.  If  the  snow  mass  is  homogeneous,  the  force  vs  sinkage  relationship 
along  the  curved  sinkage  path  should  be  very  similar  to  that  produced  by 
sinkage  along  a  vertical  or  straight  line  path.  Thus,  the  force  to  compact 
snow  along  the  curved  path  will  vary  with  distance  s  approximately  as  shown 
in  Figure  3.  The  point  labeled  ze  now  correlates  with  Sp^y  and  still 
indicates  the  point  at  which  critical  density  is  reached  and  mass  movement 
begins.  Force  at  each  point  along  the  s  axis  (previously  z  axis)  not  only 
varies  in  magnitude  but  has  a  different  direction  as  well.  In  order  to 

treat  this  situation  mathematically  it  is  required  that  an  expression  for 
force  vs  sinkage  be  available. 

Polynomial  curve  fitting  techniques  were  performed  using  existing 
force  vs  sinkage  records  generated  with  a  circular  plate-sinkage  device. 

It  was  found  that  a  fourth-order  polynomial  expression  can  be  used  to 
adequately  describe  the  force  vs  sinkage  relationship.  The  force  vs 
sinkage  expression  can  then  be  written  as 

F  =  ais*1  +  a2S^  +  a3S^  +  ai+s  +  as  (23) 

where  F  is  the  force  on  plate  and  a^ ,  a2  •••  a5  are  constants  determined 
from  a  polynomial  curve  fitting  routine. 

Substituting  from  eq  22  into  eq  23  to  get  F  in  terms  of  0 

yields 
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(24) 


F  *  ax(-4r  cos  j  +  C)4  +  a2  (-4r  cos  +  C)3 

+  a3  (-4r  cos  -|  +  C)2  +  a4(-4r  cos  j  +  C)  +  a5 

where  C  =>  4r  cos  0  a/ 2. 

This  expression  describes  the  magnitude  of  the  force  against  a  point 
on  the  tire  after  the  tire  has  rotated  6  radians  since  the  point  initially 
contacted  the  snow. 

The  direction  of  the  compaction  force  can  be  found  by  differentiating 
eq  19  to  find  the  slope  of  the  sinkage  path  at  a  given  0: 


dy  _  sinQ 
dx  1-cos 0 


(25) 


The  direction  of  the  unit  tangent  vector  t  at  any  point  on  the  sinkage  path 
can  be  expressed  as 


t 


i  + 


1-COS0  3 


1  + 


sin20  1/2 

r 


(l-cos20) 


(26) 


where  i  is  the  unit  vector  in  x  direction  (parallel  to  the  ground  surface) 
and  j  the  unit  vector  in  y  direction  (perpendicular  to  i  ). 

The  force  required  for  compaction  can  now  be  expressed  as  a  vector  (F) 
for  any  point  along  the  trochoid  describing  the  sinkage  path  from 

F  =  F  t  .  (27) 


Observing  the  tire/ snow  contact  surface,  it  is  apparent  that  each  point  on 
the  tire  surface  is  at  a  different  location  along  its  sinkage  path.  Points 
near  the  rut  bottom  have  traveled  through  most  of  their  sinkage  path.  The 
points  on  the  tire  surface  that  are  still  near  the  top  of  the  snow  have 
only  begun  travel  along  their  sinkage  path.  Thus,  each  point  on  the  tire 
is  experiencing  a  different  resisting  force.  The  sum  of  all  the  points 
along  the  tire  surface  makes  up  one  complete  sinkage  path.  If  the 
compaction  forces  are  summed  along  one  complete  sinkage  path,  the  result 
will  be  the  total  force  resisting  motion  of  the  tire.  This  can  be 
accomplished  by  integrating  force  F  over  the  distance  from  0  =  0a  (entry 
angle)  to  0  =  2ir ,  or  along  a  complete  sinkage  path.  Mathematically,  the 
total  force  resisting  wheel  motion  F-p  is 
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(28) 


_  2tt  _ 

F  =  J  F  ds 

0A 

where  ds  =  2r  sin  6/2  d0  (from  eq  22). 

Combining  eq  24,  26  and  27  into  eq  28  and  simplifying  yields 
2^ 

F T  =  /  (a^4  +  a2K3  +  a3K2  +  a4K  +  a5)  2r  sin2  -|  d6  i 

9 A  (29) 

2l1  a  q  o  fi  ft 

+  /  (axKH  +  a2KJ  +  a^YT  +  a4K  +  a5)  2r  sin  cos  d0  j 

eA 

where  K  =  -4r  cos  6/2  +  C. 


Completing  the  integration,  the  i  and  j  components  become 


-  -P* 

F  I  3 

T  i  l 


1  3  0  3  0  2  2  0  3  0 

y  cos  j  sin  j  +  — —  cos  y  sin  y 


+  (4  +  4)  cos  4  sin3  f  +  (4  +  -^)  sin3  4  (  30) 


X  X  X  X  X  X  T2t r 

-  (l6  +"8  +~)  si"  6+  +  -8+-)9je. 


T  --Pi. 

I  j  13 


2X  X 

6  0.2  5  0,3  4  0 

COS  -r  +  — r—  COS  -7T  +  — K-  COS  TT 


2  '  5 


2  2 


.  4  3  6  .  v  2 

+  — J-  cos  Y  +  X5  cos 


where  X^  =  a^  A4B 

X2  =  (4ax  A3C  4  a2A3)B 
X3  88  (6a^A^C^  4  3a2A^C  4  a3A^)B 

X^  =  (4a^AC3  4-  3a2AC^  4  2a 3 AC  4  at+A)B 

X5  ~  (axC^  4  a2C^  4  a3C^  4  a^C  4  a5)B 
A  -  -4r 

B  -  2r 

The  i  component  of  this  force  (parallel  to  the  snow-ground  surface) 
can  be  related  to  the  amount  of  resistance  that  must  be  overcome  by  the 
driving  wheels  in  order  to  advance.  The  horizontal  force  obtained  from 
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eq  30  is  acting  over  the  area  of  the  plate  used  to  generate  the 
force  vs  sinkage  curve  from  which  the  polynomial  was  acquired.  By  dividing 
this  force  by  the  area  of  the  plate  it  becomes  a  pressure  term  that  is 
independent  of  area.  The  area  of  the  tire  which  is  of  interest  for  motion 
resistance  calculations  is  the  curved  portion  of  the  tire/snow  contact 
surface.  Multiplication  of  the  horizontal  pressure  component  by  this 
contact  area  then  yields  motion  resistance  R  (for  one  wheel) ,  or 


R  = 


F  br(  A) 

T  i  A 


(32) 


MODIFIED  PLATE  TEST 

The  common  plate-sinkage  test  compacts  a  volume  of  snow  equal  to  the 
product  of  the  plate  area  and  the  snow  depth.  Following  compaction,  the 
final  volume  of  snow  is  still  of  the  same  diameter  as  the  plate  but  the 
height  has  been  reduced  by  z,  the  sinkage.  In  comparison,  a  tire  operating 
in  the  same  snow  and  at  the  same  sinkage  (same  final  volume  of  snow), 
produces  compaction  along  a  longer  sinkage  path  than  the  plate.  The 
initial  height  of  snow  compacted  by  the  tire  is  s  +  (h  -  z). 

If  an  additional  layer  of  snow  of  thickness  (s  -  z)  were  placed  on  top 
of  the  original  snow  surface  (Fig.  6)  and  a  plate  test  performed  to  a 
sinkage  of  s  (measured  from  the  new  snow  surface),  the  resulting  force  vs 
sinkage  curve  would  be  identical  to  that  produced  by  any  point  on  the  tire 
surface  if  coaxiality  of  force  and  deformation  vectors  is  assumed.  This 
curve  can  be  used  to  obtain  the  coefficients  needed  in  eq  23.  The  energy 
found  from  this  modified  plate  test  (area  under  the  force  vs  sinkage  curve) 
must  equal  the  total  amount  of  energy  required  to  advance  the  tire.  This 
can  also  be  called  the  total  energy  of  compaction  and  includes  both  the 
vertical  and  horizontal  sinkage  components. 

During  compaction  with  a  tire,  the  snow  is  vertically  translated  for 
a  distance  equal  to  z,  the  compaction  which  takes  place  during  a  standard 


s 

L  Figure  6.  Comparison  of  sinkage  paths 
_  during  compaction  by  a  plate  and  a  wheel. 
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plate-sinkage  test.  The  energy  found  from  this  test  can  be  equated  to  the 
vertical  component  of  compaction. 

If  the  total  energy  of  compaction  along  the  curved  sinkage  path  (from 
the  modified  plate  test)  and  the  vertical  component  (from  a  standard  plate 
test)  are  known,  the  horizontal  component  (the  energy  required  to  translate 
the  snow  horizontally)  can  be  computed.  This  energy  can  also  be  considered 
as  the  sum  of  the  horizontal  force  components  acting  on  the  tire/ snow 
contact  surface.  Thus,  vehicle  motion  resistance  for  one  wheel  (per  unit 
of  distance  traveled)  can  be  found  from  plate  tests  in  the  same  snow  from 


R  = 


br(2tr  -  0A) 


A 

P 


(33) 


where  is  the  energy  (or  sum  of  forces)  determined  from  a  modified 
plate  test,  and  Eg  the  energy  (or  sum  of  forces)  determined  from  a 
standard  plate  test. 

Although  this  discussion  has  chosen  to  represent  the  force  vs  sinkage 
relationship  for  snow  by  a  fourth-order  polynomial ,  the  same  process  of 
breaking  down  snow  compaction  into  vertical  and  horizontal  components  can 
be  accomplished  with  other  mathematical  relationships.  In  some  instances, 
logarithmic,  hyperbolic  or  exponential  curves  may  best  fit  the  data 
produced  by  a  modified  plate  test.  By  redefining  F  (as  a  function  of 
sinkage  s)  in  eq  23,  vertical  and  horizontal  force  components  can  be 
computed  from  eq  27  and  28. 

Conditions  of  positive  and  negative  slip  can  also  be  taken  into 
account  by  not  equating  a  and  r  in  eq  19.  This  leads  to  the  following 
expressions  for  distance  along  the  curved  sinkage  path  and  the  unit  tangent 
vector: 


s 


2*tt 


(a^  +  r^  -  2  ar  cos  0)*^  d0 


(34) 


t 


.  r  sm  0 

l  +  - -  j 

a  -r  cos  0 


r  sin  0 
a  -r  cos 


T U2~‘ 
■)  1 


(35) 


Substituting  these  expressions  into  eq  23  and  27,  and  completing  the 
integration  in  eq  28,  yields  values  for  the  vertical  and  horizontal 
components  of  compaction. 
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APPLICATION 


A  simple  test  was  performed  to  get  an  idea  of  the  shape  of  the 
deformation  area  ahead  of  a  tire  traveling  in  snow.  Several  closely  spaced 
vertical  holes  were  punched  in  a  virgin  snow  underlain  by  concrete. 

The  holes  were  each  dusted  with  colored  chalk  powder  to  coat  the  sides  of 
the  hole.  Approximately  half  of  the  holes  were  then  overrun  by  a  tire.  A 
vertical  section  sliced  through  the  holes  revealed  the  deformed  area  (Fig. 
7).  It  is  readily  apparent  that  a  great  deal  of  horizontal  displacement 
has  taken  place.  It  can  also  be  seen  that  the  basic  outline  of  the 
deformation  bulb  is  as  depicted  in  Figure  4b. 

The  opportunity  to  validate  the  concept  of  the  assumed  coaxiality  of 
stress  and  deformation  vectors  has  not  yet  presented  itself.  It  is 
concluded,  from  theoretical  concepts  (Richmond  and  Blaisdell  1980)  and 
observed  tests  (both  plates  and  tires  in  snow),  that  this  concept  is 
sound.  Field  and  laboratory  tests  must  be  performed,  however,  to  assess 
the  accuracy  of  dividing  the  energy  of  compaction  into  vertical  and 
horizontal  components  as  shown  in  eq  30  and  31.  These  tests  must  also  be 
used  to  determine  the  accuracy  of  eq  33  and  the  modified  plate  test 
concept* 

In  its  present  form,  vehicle  motion  resistance  can  be  predicted  for  a 
given  snow  by  performing  the  modified  plate  test  described  above  in  the 
snow  and  determining  the  coefficients  of  the  best  fit  fourth-order 


Figure  7.  Chalk  dust  holes  depicting  deformation 
area  and  compaction. 


16 


polynomial.  These  coefficients  and  the  vehicle  sinkage  can  then  be  used 
with  eq  30  and  32  to  predict  R. 

An  alternate  approach  (which  should  yield  approximately  the  same 
value)  is  to  run  both  a  modified  and  standard  plate  test  in  snow  and  use 
their  force  vs  sinkage  curves  with  eq  33  to  predict  R. 
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